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Abstract 

Let {g n amb)n,m& be a Gabor frame for L 2 (M) for given window g. 
We show that the window h° = S~~?g that generates the canonically 
associated tight Gabor frame minimizes \\g — h\\ among all windows 
h generating a normalized tight Gabor frame. We present and prove 
versions of this result in the time domain, the frequency domain, the 
time-frequency domain, and the Zak transform domain, where in each 
domain the canonical h° is expressed using functional calculus for 
Gabor frame operators. Furthermore, we derive a Wiener-Levy type 
theorem for rationally oversampled Gabor frames. Finally, a Newton- 
type method for a fast numerical calculation of h° is presented. We 
analyze the convergence behavior of this method and demonstrate the 
efficiency of the proposed algorithm by some numerical examples. 
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1 Introduction 



A Gabor system consists of functions of the form 

Qna^it) = e 2 ™ mbt g(t - no), n, m e Z, (1) 
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for some window g G -L 2 (IR) and a, b G M. The parameters a and b are the 
time- and frequency translation parameters. 

We say that the triple (g,a,b) generates a Gabor frame for 2v 2 (R) for 
given a, b if there exist constants (frame bounds) A,B>0 such that 

A||/|| 2 < \(f,9na, mb )\ 2 <B\\f\\ 2 , (2) 

for any / G L 2 (R). 

The Gabor frame operator S is given by 

Sf = (f'9na,mb)g na,mbi f e L 2 (R), (3) 

and satisfies 

AI < S < BI, (4) 

where I is the identity operator of Zv 2 (IR). 

If {9na,mb)n,mez is a Gabor frame for 2v 2 (M) then any / in X 2 (1R) can be 
written in the form 

/ ^ ^ (/) 1 na,rnb) 9na,mb ^ ^ (/; 9na,mb)'~J na,mbi (5) 

where the canonical dual frame (lnamb)n,m€Z is given by 

< a ,mb = e 2mmb V(t-na), n,meZ (6) 

with 7° = S~ 1 g. 

If A = B the frame is called tight. In this case S = -jl, whence 7 = j^g. 
We say that a frame is normalized tight if in addition A = B = 1. Tight 
Gabor frames play an important role in signal processing and communica- 
tions. They appear for instance in the construction of paraunitary modulated 
filter banks |], |J as well as in the construction of orthogonal frequency di- 
vision multiplex (OFDM) systems || |32| in wireless communications. Tight 



Gabor frames are also useful tools for the analysis of pseudodifferential op- 
erators 127 . 



Given a frame (g n a,mb)n,mez with frame operator S a standard technique 
to construct a tight frame is the following. Set 

h° = S-ig, (7) 
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then it is easy to see that {h^ n a,mb)n,m& is a tight frame for 2v 2 (IR). We will 
refer to h° as canonical tight window and to (/ina,mb)n,mez as canonical tight 
frame. 

It is known that the canonical dual 7 minimizes 



9 7 



It 



(8) 



among all dual windows 7, cf. [|L9 |. This minimality property of 7 can be 



expressed more precisely in the time-domain as: for a.e. t and all a > 0, f3 > 
(possibly depending on t) we have that 

00 

\^9(t - na) - 07(t -na)\ 2 (9) 

n=— oo 

is minimal among all dual 7 for 7 = 7°. Such minimality results also hold in 
the frequency domain and in the Zak transform domain. 

We will show in Section ^| that a result similar to (g) holds for tight Gabor 
frames. More precisely, we will demonstrate that h° solves 

argmin||g — h\\ (10) 

heL 2 (R) 

among all h such that (h, a, b) is a normalized tight frame. 

This result is also important from a practical viewpoint, where one often 
selects a "nice" g and constructs the tight frame by using (0) with the goal 
in mind that the resulting tight frame window should have similar properties 
as g. Thus it is gratifying to know that h° is closest to g (in the i 2 -sense) 
among all possible normalized tight frames (h,a,b). For conditions under 
which h° inherits the decay properties of g see |35], |31], |5), [!2| . 

It is sufficient to restrict our attention to normalized tight frames due to 
following reason. Assuming that (ITO) holds the reader may easily convince 



herself /himself that the window h° that minimizes \\g — h\\ among all h such 
that (h, a, b) is a tight frame is given by 

(h°,h°) [ ' 

For the special case that (g,a = 1,6 = 1) constitutes a frame it is not 



difficult to see that (10) holds. In this case we can consider the problem in 
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the Zak transform domain with A = p = q = 1 (see Subsection |1.1| (|2~3|)-(2T)) 
for the definition of the Zak transform and the Zibulski-Zeevi matrices that 
are implicitly used in the next steps). Thus we have for any / such that 
(/, 1, 1) has a finite frame upper bound that 



(Z 1 Sf)(t,u) = \{Z 1 g){t^)\\Z 1 f){t,iy), a.e.t,zA 



By functional calculus (cf. Subsection 

(Z 1 h°)(t,u) = (Z 1 S-$g)(t,i') = 

(Z x g)(t,v 



II) we then have that 



Z l9 )(t,v)\ 2 ) \Z l9 )(t,u) 



a.e. t, v. 



\(Zig)(t,v) 

Moreover, normalized tightness of (h, 1,1) is equivalently expressed as 

\(Zih)(t, u)\ — 1, a.e. t, v. 
Now when we have a complex number z then 





z 




\z\ - 1 


min \z — uj\ = 








\z\ 







Thus we have that for a.e. t, v and all normalized tight (h, 1, 1) 



(Z 1 g)(t,v)-(Z 1 h°){t,u) 



< 



[Z l9 ){t, u)-{Z 1 h){t,u) 



(12) 

(13) 
(14) 

(15) 
(16) 

(17) 



Integrating the latter inequality over [0, l) 2 we get using unitarity of Z x 
-> L 2 ([0, l) 2 ) that \\g - h°\\ 2 < \\g - h\\ 2 . 



The rest of the paper is organized as follows. Subsection |1.1| contains a 
brief review of various representations of the Gabor frame operator in differ- 
ent domains. We will make use of these representations in Section |2| when 
we prove the minimization result (0). We also discuss some refinements and 
extensions of this result. In Section [|we present a Wiener-Levy type theorem 
for Gabor frames. Finally, in Section ^ we describe a Newton-type method 
for computing h° and analyze the convergence properties of this method.Q 

It is obvious that some or many of our results admit generalization to 
more general lattices in IR 2 as sets of time-frequency points to which the 
windows are shifted and to higher dimensions, but we have not pursued this 
point here. 

1 We were kindly informed by Hans G. Feichtinger that he discovered the algorithm of 
Section ^ independently of us, as early as 1995 and used it in numerical experiments to 
compute h° with satisfactory results. 
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1.1 Representations of the Gabor frame operator 

We briefly review the representations of the Gabor frame operator in the 
time domain, in the frequency domain, in the time-frequency domain, and, 
for rational oversampling where ab G Q, in the Zak transform domain. We 
also present a lemma concerning the approximation of the inverse square root 
of an infinite-dimensional matrix by finite-dimensional inverse square roots. 

Assume that (g, a, b) has a finite frame upper bound B, and define the 
Ron-Shen matrices 

MJt) := (g(t-na- k/b)) , a.e.t, (18) 

V / fcez.nez 

H g (v) := [g(u — mb — l/a) ) , a.e. v. (19) 

Now there holds that M g (t) and H g (v) are for a.e. t and a.e. u bounded linear 
operators of £ 2 (7j) (with operator norms < (bB)^ and < (0.8)2). Now when 
(/, a, b) also has a finite frame upper bound, then there holds 

M Sf (t) = ^M g (t)M* g (t)M f (t), a.e.t, (20) 

Hsfiy) = ^H g (u)H* g (u)H f (u), a.e. v. (21) 



For the case (|T9[) , (PH) this can be found in |2TJ, proof of Theorem 1.2.6 and 



Proposition 1.2.3, in the more general context of shift-invariant systems, and, 



in somewhat different form, in [|28[ or in [29]. 

Next, still assuming that (g, a, b) has a finite frame upper bound B, we de- 
fine the analysis operator with respect to the adjoint lattice {(k/b, l/a)\k, I G 
Z}by 

m G L 2 (R) -> U g m = ({m, g k/b ,i /a }) • (22) 

Now there holds that U g is a bounded linear operator (operator norm < 
(abB)^). When (f,a,b) also has a finite upper frame bound, then there 
holds 

U Sf = ^- b U 9 U* g U f . (23) 



These results can be found in ||19|| , Propositions 2.7, 3.1, and also in |29 
and (in somewhat different form). Expression ( |23| ) is sometimes called 
the Janssen representation of the frame operator. 
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Finally still assuming that (g, a, b) has a finite frame upper bound B, 
but now with ab = p/q G Q and 1 < p < q, gcd(p,q) = 1, we define the 
Zibulski-Zeevi matrices as 

&(t,v)= p -h((Z ig )(t (24) 

\ b IP ' k=0,...,p-l;l=0,...,q-l 

W{t, v) = p~* ((Z a g)(t -k±v- -)) (25) 

V PI ' k=0,...,p-l;l=0,...,q-l 

for a.e. t, v G M. Here we define for A > the Zak transform Z\f of an 
/ G L 2 (R) as 

oo 

(Z x f)(t, u) — A^ ]T /(A(t - A;))e 2 ^, a.e. t,i/6i (26) 

fc=— oo 

Now there holds that v) and ^ 9 (t, v) are for a.e. t, v bounded linear 
mappings of C p into C q (matrix norms < B?). When (f,a,b) also has a 
finite frame upper bound, then there holds 

§ s f(t,v) = $ 9 (t,z/)(V(t,z/)) V(t,z/), a.e.t,z/, (27) 
\S? s f(t,p) = ^ 9 (t,z/)(V(t,z/)) V(t,z/), a.e.t,z/. (28) 

These results can be found in [^(J and in pj|j . 

The frame operator S is represented through the operators 

S(t) ■= ±M s (t)M'(t), S(v) := iff s Mfl», (29) 

a==s<w < 30) 

and 

5 x (t,i/) :=$ 9 (t,z/)(V(^))*, 5 2 (t,i/) !/)(*'(*,!/))* (31) 

in the following; sense. When A > 0, B < oo, g G £ 2 (M), then 

(g, a, 6) is a Gabor frame with frame bounds A, B (32) 
-<=>- A/ < S(t) < BI,&.e.t 
^ AI < S(y) < BI,a.e.u 
AI <S<BI 

AI < Si(t, v) < BI, a.e. t, v 
AI < S 2 (t, v) < BI, a.e. t, v. 
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Here the J's denote identity operators of £ 2 {Z), £ 2 {Z), £ 2 (Z x Z), C p , and C p , 
respectively. See the references just given. 

An important point for the developments in this paper is that the relations 
(|20|), (pl|), (p3|), (p7|), (p8|) can be extended as follows. Assume that (g,a,b) 
is a frame with frame bounds A > 0, B < oo, and let 9? be a function analytic 
in an open neighbourhood of [A, B]. When (/, a, b) has a finite frame upper 
bound, there holds 

M vm (t) = (p(S(t))M f (t), a.e. t, (33) 

H v{8)f (u)=ip(§(u))H f (u) 1 a.e.1/, (34) 

^(5)/ = ¥>(£)£//, (35) 

§vWf(t,u) =(p(S x (t,u))$ f (t,u), &.e.t,v (36) 

=p(S 2 (i, 1/), a.e.t,!/. (37) 

Indeed, the aforementioned relations can be shown to hold by iteration when 
if is a polynomial, and then a standard argument in functional calculus not- 
ing ( |3^ ) then yields the result for general (p analytic in an open neighbour- 
hood of [A, B]. In the important special case that (p(s) = s a , one can do this 
extension explicitly by using 

s ^(^ {1 - v) y = (E±Ayy(-) hvr , (38) 



n=0 



n 



where v = l-2s(B + A)' 1 has modulus < (B- A){B + A)~ l when s G [A,B\. 

A further result that we shall use is that for any frame (g,a,b) and any 
normalized tight frame (h, a, b) we have 

(g,S- 1 g) = \\h\\ 2 = ab. (39) 

See H H, gjj|. 

In the sequel we will also make use of the following result. 

Lemma 1.1 Let T = {Tki)k,i& be a hermitian positive definite biinfinite 

_l _i 

matrix with c\I < T < c%l . Define Tjv = (Tki)\k\,\i\<N- Then T N 2 — > T 2 
strongly for N — > 00 . 
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Proof: Without loss of generality assume c?, < 1, otherwise we can always 
set V = —T. Note that 



C2 



r-* = £>(/-T) fc , (40) 

fc=0 

where c fc = (7) and T" 2 " = ^Zo - T N f . Define = £f =0 c fc( J - 
T) fc and rjjf ) = Ef=o c fc( J JV - ^v) fc - Observe that 

C \I — \\Tn\\ — C2I uniformly in N (41) 

and 

T N ^ T strongly for N -> 00. (42) 

For fixed iV there holds Tjv — > T^ 2 as X — > 00. Furthermore, for each K 
we have that for N 00. This together with properties (41) 

and (pE2|) implies that T^ 2 — > strongly for N — > 00. □ 

2 Canonical tight Gabor frames and the 
Fan- Hoffman inequality 

Theorem 2.1 Let (g,a,b) be a frame and let h° = S~^g be the canonically 
associated tight frame generating window. Then for any h E L 2 (K) for which 
(h, a, b) is a normalized tight Gabor frame there holds 

< \\g-h\\ < \\g + h°\\. (43) 



In Subsections |2.1] , [2.2| , and |2.3| we will present different proofs of this 
inequality, making use of different representations of the Gabor frame oper- 
ator. 



Inequality ( |4*3| ) can be considered as a Fan-Hoffman type inequality JTT ] . 
The Fan-Hoffman inequality can be formulated as follows. Let the matrix 
M G C mxn , m > n have the polar decomposition M = UH, where U*U = I n 
and H is hermitian positive semidefinite. Then for any mxn matrix W with 
W*W = L 



\M - U\\ < \\M - W\\ < \\M + U\\ (44) 
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for any unitarily invariant norm. See Jl ]J for a proof for M G C raxn and [T7 
for a discussion of the case M G C mxn . 

An immediate consequence of the Fan- Hoffman inequality is the following 
result. Consider a frame for C n consisting of m vectors (fk, k — 1, . . . , m with 
frame operator S. Then the tight frame (/0fc)fcLi that minimizes 



£ 

k=l 



Vk-^kW ( 45 ) 



among all tight frames (V'fc)fcLi hi C n is given by 



, i 



trace(52j i 
= ^ — *(p k . (46) 

This can be easily seen by considering the vectors (fk as rows of an m x n 



matrix M and using ( pE4"D with the Frobenius norm, see also Section 1.5 in p3 

The case where n = m and {fk)k=i ^ s an orthonormal basis has been in- 
vestigated by Lowdin in the context of quantum chemistry [16| . The solution 



- known under the name Lowdin orthogonalization - is of course a special 
case of fl4"S|). 

Both, the Fan-Hoffman inequality and the Lowdin orthogonalization, 
have been extended to infinite dimensions, see e.g. Chapter VI in [15|, or Jl|, 



13 . These extensions usually require that the involved operators satisfy some 



Hilbert-Schmidt-type condition in order that 

and the operator version of (fPj) are finite. A different notion of closeness of 



frames, based on work of Paley and Wiener [ 26| , has been used in . 

For Gabor frames the expression J^mnez \\9na,mb ~ h na)mb \\ 2 is of course 
not finite for g ^ h. But since \\g n a,mb — h na ^ m b\\ = \\g — h\\ for all n,m 
we can restrict ourselves to the minimization of ||g — h\\. The minimization 
result for tight Gabor frames presented in this paper can also be considered 
directly in the context of Lowdin orthogonalization, due to the following 
reason. The results in Subsection |1.1| imply that if (h, a, b) is a tight frame 
for L 2 (R) then (h,l/b,l/a) is an orthogonal basis for its closed linear span. 
Thus the statement that h° minimizes ||g — h\\ among all normalized tight 
frames (h, a, b) is equivalent to saying that h° minimizes \\g — h\\ among all 
h such that (h, 1/6, 1/a) is an orthonormal basis for the closed linear span of 

i,9na,mb)n,m- 
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2.1 The problem in the time domain and in the fre- 
quency domain 

2.1.1 The problem in the time domain 



For the first proof of Theorem |2.1| we consider the Gabor frame operator in 
the time domain using Ron-Shen matrices. 
We consider matrices 

MJt) = ( g(t - na - r)) a - e -^ ( 48 ) 

V b J kez,nez 

with row index k and column index n. We have by functional calculus, 
see (H) 

M h o(t) = (-M g (i)M;(i)J a.e.t. (49) 

The condition of normalized tightness can be expressed equivalently as, 
see (|32l) 

M h (t)M* h (t) = bl a.e.t. (50) 
Note also that by normalized tightness \\h\\ 2 = ab (see ([39)). Thus 

\\h - g\\ 2 = ab+ \\g\\ 2 -2Re{h,g), (51) 

and 

\\h + g\\ 2 = ab + \\g\\ 2 + 2 Re{h, g), (52) 

and so we are to maximize Ke(h,g) over all h satisfying floTI). 
Let K = 1,2,... and consider the sections 

M g>K (t)=(g(t-mi-j)) . (53) 

V b J\k\<K,neZ 

M g ^K consists of (2K + 1) rows and a biinfinite number of columns. Note 
that 

M g Mt)Kx® = ((^W^W)^!*!^,^- (54) 
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Thus we have that for any allowed h there holds 

M g>K {t)M* g>K (t) = bI 2 K + i. (55) 
We are going to maximize, for fixed t £ R, the quantity 

trace (M K M* tK (t)), (56) 
where Mk = (Mk n )\k\<K,nez is a matrix satisfying 

M K M* K = b (57) 



Noting that the matrix in (J54]) is positive definite (as a (2K + 1) x (2K + 1) 
section of the positive definite operator M g {t)M*{t) since (g,a,b) is a frame 
for L 2 (R)) we do a singular value decomposition of M g ^{t) according to 

2K 
3=0 



where 

^ T £ € 2 ( 

IT r^A+l 



^ = fen)„ 6Z e ^(Z) 



U,- = [Uj-K, • • • , £ 

are orthonormal, <r| is the j-th eigenvalue of M 9i K(t)M* K (t) and «j is the 
corresponding eigenvector. Thus 

(M g , K {t)Ml K {t)) Uj = a)ui (59) 

and 

v 3 = ^-Ml K (t)u r (60) 
a 3 

Now we compute 



2A' 



trace (M K M* g K {t)) \ = | trace (M K ^ a j v j u f) I 

3=0 

2K 2K 

trace ( UjMkVjuf) | = | trace ( <jjuf MkVj) | 

3=0 3=0 
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2K 2K 

^2 °o ( m kVj ,Uj)\ < y^QjH M K V 

3=0 3=0 



j ii ii M i I 



2K i 

< b 1 ? a 3 = 6 * trace ( M 9,x(t)M*^(t)) 2 . (61) 



3=0 



Taking 



M K {t) = (-M g>K {t)M; iK (t)} *M g>K {t) (62) 
we also see that the maximum value 

i 

65 trace ^M gtK (t)M* gjK (tyj 3 (63) 



is assumed by M K {t) in (j62|) . 

Next we observe that for an allowed h we have 

K oo , . 

tm C e[M KK {t)M; K {t)] = h ( f ~ na ~ - no - -). (64) 

k=—K n=—oo 

Taking the average over and integrating over t G [0, a), we thus see that 
(h, 9) — J ££=-* W - na - !)<?*(* - na - f ) dt 

° a (65) 

= 2W+1 1 trace[M h , K (t)Ml K (t)]dt. 


The right-hand side of fl65D has modulus < what one would obtain by re- 



placing M h ^ K (t) by M K {t) in (|65|). Due to Lemma O] we have that 



lim M K {t) = lim (\M 9tK {t)M*{t)) * M 9tK {t) = M h o(t) (66) 

K— >oo K^oo\0 a ' / 



in the strong operator topology This implies that 



a 

I ^i^^KK^dt = (h°,g), (67) 
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and the proof of inequality ( fj3|) is complete. 



Actually we have shown more than fl43|). Consider the a-periodic L°°- 
functions 

/(*)= Ht-na)g*(t-na) = (h,g) a (t), (68) 

n=— oo 

where we have used the bracket product notation (•,•)<, as in [0]. That 
/ G follows from the frame bound conditions, together with the con- 
siderations in Subsection |ET| . Consider the case that ab = p/q with integer 
p, q, 1 < p < q, gcd(p, q) = 1. Then for a.e. t there holds (see the right-hand 
side of 



2 At £ /(t - (69) 



k=-K 
K 



y f{t- q -a) 

1 ; — 



fc=-iC 

[fc=W] 2iT + 1 2 S f(t-lqa- 

r=0 l;-K<lp+r<K y 



I P_1 g 
— — r / . #i-K<i P +r<K ■ fit — r-a) 



[using the a-periodicity of /] 27^ -|- 1 ' ^ ' p 

1 + 1 ^ 



2K+1 p 
1 

±$>,0>.(*-r/&) 



r=0 

p-1 



^ r=0 

Thus there holds, more precisely, that for a.e. t 

p-i 

P 



-J2(h,g) a (t-r/b) (70) 



r=0 



is maximized among all normalized tight (h, a, b) by h°. And then, of course, 
integrating over t G [0, 1/6) gives the result. 
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In case that ab ^ Q, it can be shown that 

— i— f(--k/b)^(h,g) (71) 

k=-K 

in the mean square sense when K— > oo (here / is as in fl68[ )). Hence in this 
case we do not get a more precise result of the type given above for rational 
ab. 



2.1.2 The problem in the frequency domain 

We now present a proof of Theorem ^J] where we consider the Gabor frame 
operator in the frequency domain, using frequency domain Ron-Shen matri- 
ces. We thus consider 

H„(u) = [g(u — mb — l/a)) , a.e. v, (72) 

and we have by functional calculus 

H h o(v) = (-H g {v)H;{v)Y~ 2 H g {v), a.e.i/, (73) 

while the condition of normalized tightness can be expressed equivalently as 

H h (u)H* h (u) = aI, a.e. v. (74) 

We therefore find ourselves in a situation that is similar to what we had in 
Subsection |2.1.1| in all respects, except that the Ron-Shen matrices involve 
now the Fourier transform of the windows, rather than the windows them- 
selves. Accordingly, we know that \\h — g\\ 2 is minimal among all h such that 
(h, a, b) is normalized tight for h = h°. By Parseval's theorem we have that 
\\h — g\\ 2 = \\h — g\\ 2 , and the proof is complete. We also get a frequency 
domain sharpening of the result for ab G Q in a similar manner as the one 
obtained at the end of Subsection [2.1.1| . This latter sharpening reads that 
for a.e. v 

-T(h,gW~s/a) (75) 
is maximal among all h such that (h, a, b) is normalized tight for h = h°. 
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The frequency domain result can be generalized to some extent to shift- 
invariant systems 

(gn,m)n,meZ = \gm(- - na)j , (76) 

where g m E L 2 (JBL) and a > 0. Following the presentation of the Ron-Shen 
theory |28[] given in [^TJ, Subsection 1.2, we now consider the frame operator 



Sf = y 52(f,g n , m )g n , m , feL 2 (R), (77) 

n,m 

and the Ron-Shen matrices 

H g (u) = (g m (u -l/a)) , a.e. v. (78) 

There holds for A > 0, B < oo that 

AI < S < BI <=> AI < -H g {v)H*(v) < BI, a.e. v. (79) 

The canonical dual system (7° m ) n ,mez an d t ne canonically associated tight 
system {h\ m ) n , m&L are given by 

7 ^(. -na)) , (h° m (--na)) (80) 

with 7^ = S~ l g m ,h? m = S~^g rn . Here we note that 5 in ( [771 ) commutes 
with all shift operators /—►/(• — na),n G Z. The conditions of duality 
of (7nm)n,m6Z and normalized tightness of (h nTn ) n)m ez can be equivalently 
expressed as 

H g {v)H*(v) = al, H h (u)H*(u) = al. (81) 

We have a functional calculus in terms of the matrices a~ l Hg(v)H*(y) by 
the formula 

H v{S)g {v)= v (±H g {v)H;{v))lI a {v), a.e. i/, (82) 

where 99 is a function analytic in an open neighbourhood of [A, B}. We have 
set here 

( fi{S)g = ({(p{S)g m )(- - na)) , a.e. v. (83) 
15 



In particular there holds 

H h o(u) = (~H g {v)H* g {v)) * H g (u), a.e.z/. (84) 

Now the developments as given in Subsection [2.1.1| can be mimicked to a 
large extent. Considering truncated matrices H gK as in (|53|), and denoting 
the Frobenius norm by \\-\\f, the minimization of (for a.e. v) 

oo 

\\E K -E 9)K (y)f F =Y, E \h m -g{y-lla)\ 2 (85) 

\l\<Km=-co 

over all (2K + 1) x oo matrices Hk = (hi m )\i\<K,mez satisfying 

H K H* K = aI 2K+1 (86) 
yields that the minimizing Hk is given by 

H K = (±H g , K (u)Hl K (v)y*H g , K (v). (87) 



And, again by Lemma |1 . 1| , this Hk — > H^o as K— > oo in the strong operator 
topology. 

2.2 The problem in the time-frequency domain 



We next present a proof of Theorem [2.1| where we consider the Gabor frame 
operator in the time-frequency domain, using analysis operators U g as defined 
in (p2|) based on the adjoint lattice {(k/b, l/a)\k, I E Z}. 
Note that U g U* has the matrix representation 

U M=(<9*M/a,9k M *))^j, & - (88) 
Normalized tightness is equivalently expressed as 

U h U* h = abI. (89) 
The canonical tight frame generating h° is given by 



1 

ab 



u h° = ( —u g u; ) u g (90) 
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by functional calculus. We also note that 

u h u: - 



Now consider the sections 
U g jc, L : m G L 2 (R) — > 

with adjoint s 



(9k'/b,l'/a, hk/b,l/a)j 

(m, gk/b,i/a) 



fc.JeZ;fc' i'ez 



e C 2im x C 2L+1 



\k\<K,\l\<L 



u 



g,K,L 



: c G C 2A ' +1 x C 2L+1 - C ^/» 



l/a 



\k\<K 

\1\<L 



and 



U g ,K,LU* KL - ({Ug tKjL U* iK)L ) 



\<K,\l\<L;\k'\<K,\l'\<L 

This we also do with allowed /i's, see (189D, and we thus see that 



Uh,K,LU^ K ,L - abI(2K+l)(2L+l)- 

Now consider the maximization of 



over all 



U. 



K,L 



trace 



Uki-k'v 



Uk,lU*k,l 



\<K,\l\<L;\k'\<K,\l'\<L 

such that Ur,lU* kl = abI(2K+i)(2L+i)- This yields the maximum value 



trace 



Ob)*[Ug,K,LU* g ,K,L 



assumed by 



Note that from 



U 



K,L 



1 



-U nK rU*, 



'U n 



(g,h) 



ab ~g,K,LVg tKtL ) Ug )K>L . 
1 



(2K+ \){2L + l 



trace 



UhK,lU* 



g,K,L 



By Lemma |TTT] Uk,l in ( p9|) converges strongly to Uho : K,L, hence 

1 



(9,h")= lim 



k,l^oo(2K + 1)(2L + 1 
Therefore, |(<7, < (g,h°). 



■ trace 



Uk,lU*,k,l 



(91) 
(92) 
(93) 

(94) 

(95) 
(96) 
(97) 

(98) 
(99) 
(100) 

(101) 
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2.3 The problem in the Zak transform domain 



We proceed by presenting a proof of Theorem |2.1| by considering the Gabor 
frame operator in the Zak transform domain, using Zibulski-Zeevi matrices. 
The proof in this section is the simplest among all presented proofs, however 
it only applies for rational oversampling. 

Assume ab — p/q, 1 < p < q,gcd(p,q) = 1. By J2TJ, Sec. 1.5, (1.5.11) we 
have 

q ~ lp ' 1 p-l q-l 

J J / n 7 n 







k=0 1=0 



I P 

p J J trace [$ h (t, v) ($ 9 (t, u))*] dt dv. 







Here $^(t, v) is the Zibulski-Zeevi matrix as defined in fl24|). The condition 
of normalized tightness is expressed equivalently as 

<£> h (t,v)(® h (t,v))* = I p &.e.t,v. (102) 

The maximum of 

trace [$($ 9 (t, z/))*] (103) 
over all p x q matrices $ with <3>$* = I p is equal to 

trace [(V(i, u) ($ 9 (t, z/))*) "1 (104) 

and assumed by 

_ i 

$(t,u) := (® 9 (t,is)($ 9 (t,is))*) ~ 2 1> 9 {t,v). (105) 



The right-hand side of (|105|) equals $ho(t, v) by functional calculus for Zibulski- 
Zeevi matrices. 

Hence we see that h° indeed maximizes (h,g), and that it does so in a 
pointwise manner in terms of Zibulski-Zeevi matrices. 

The proof using the matrices as defined in ( p5|) is almost identical to 
the one for the matrices $ 9 and is therefore left to the reader. 
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3 A Wiener-Levy theorem for Gabor frames 



In [31], [5|, |32| it has been investigated under which conditions the canonical 



tight window inherits the decay properties of the window g. Loosely spoken 
h° inherits the decay properties of g whenever this is true for 7 . We give 
a precise mathematical formulation of this observation. We start with the 
following result from the theory of Banach algebras |T| . 



Proposition 3.1 Let B be a Banach algebra, and x an element of B with 
real-valued spectrum a(x) > 0. Then there exists in B a unique square root 
of x with real-valued positive spectrum. 

This implies the following: Assume that the properties of g imply that 
S belongs to some Banach algebra B. Whenever S~ x G B, it follows that 
S~* G B. 

Actually for rational values of ab a more general result holds. 

Theorem 3.2 Assume that ab = p/q,l < p, q G Z,gcd(p,g) = 1. Assume 
that (g,a,b) is a frame, where g satisfies Tolimieri and Orr's condition A 

(cfWM) 

^2\(9,9k/b,l/a)\ < OO. (106) 

k,l 

Let W be a time-frequency shift-invariant Banach space of windows w such 
that (w, a, b) has a finite frame upper bound when w G W. Finally, let (p be 
analytic in an open neighbourhood of [A, B], where A > 0, B < 00 are the 
frame bounds of (g, a, b). Then ip(S)w G W for any w G W. 



For the proof of Theorem |3.2| we first need a matrix- version of the classical 
Wiener-Levy theorem. We introduce the following 

Definition 3.3 A set U of functions f(t) = Ylkez a k{f) e2mkt ^ s sa ^ ^° have 
an absolutely uniformly convergent Fourier series if for every e > there is 
a K > such that 

^Klf)|<£, for all feU. (107) 

\k\>K 
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Proposition 3.4 Let H(t) be a periodic hermitian matrix-valued function 



k 



with absolutely convergent Fourier series, i.e., H(t) = ^2 ke %Hke 2 t ,H l 

((#fc)m,n)m,„=l and 

\\H\\ T := ^i™ m ,n|<oo. (108) 



N 



m,n=l kalh 



Let A, B G K, A < B be such that all eigenvalues of all H(t) are in [A, B\. If 
the function (p(z) is analytic in an open neighbourhood of [A, B] then if(H(t)) 
also has an absolutely convergent Fourier series. 



Proof: We have the Dunford representation Jit 

<p(H{t)) = ^~j (zl - Hit))' 1 dz, (109) 
r 

where T is a closed contour containing [A, B] in its interior. We write 

(zl - H(t)) ^ = det (zl -Hit))' 1 adj (zl -H(t)). (1 10) 

The entries of zl — H(t) have Fourier series absolutely uniformly convergent 
in z G T. Hence it follows easily that adj (zl — H(t)) and det [zl — H(t)) , as 
sums and products of the entries of zl — H(t), have Fourier series that are 
absolutely uniformly convergent in z G T. Furthermore, there are constants 
C\ > 0, C 2 < oo such that 

C x < | det (zl - Hit)) | < C 2 , zeT,teR. (Ill) 

Thus by the uniform Wiener l//-theorem we have that det (zl — Hit)) 1 
has a Fourier series that is absolutely uniformly convergent in z G T. This 
implies that — H)~ 1] \f is uniformly bounded in z G T, by C, say. Hence 

\\<p(H)\\r< j \y{z)\\dz\ < oo, (112) 
r 

where the boundedness of the integral in ( |112j ) follows from the continuity of 
if. The proof is complete. □ 
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The knowledgeable reader will have no difficulties to extend the argument 
concerning the uniform Wiener 1 / /-theorem given in || to the case that we 
are considering functions of two variables with absolutely uniformly conver- 
gent Fourier series. With the two-dimensional version of the uniform Wiener 
l//-theorem at hand it is then easy to extend Proposition |3.4j to bivariate 
hermitian matrix- valued functions H(t, v). 

Proof of Theorem |3.2| : Let w G W. By functional calculus we have 



kl 



0,0;k,l 



Wk/b,l/ai 



(113) 



with convergence of the right-hand series at least in weak X 2 -sense. It is 
therefore sufficient to show that 



kl 



0,0;fc,Z 



< OO 



(114) 



since Wk/b,i/a G W for all k, I G Z. 

To show ( 114j) we take w = g in ( |113| ). By biorthogonality we have 



From Prop. 1.1 in |2(J we have for any /, h G L 2 (R) 

^2\{f,h k /b,i/a)\ < oo 

k,l 

$^(£, v) (& h (t, u)Y has an abs. conv. Fourier series. 

Now take 

/ = <p{S)g> h = 1 ° = S-'g. 
By functional calculus in the Zak transform domain, we have 

= <J$ 9 (t, v){$ 9 {t, u))' 



(115) 



;ii6) 



;ii7) 



118) 
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By the bivariate version of Proposition O we have that the right-hand side 
of (|118[) has an absolutely convergent Fourier series. Then by Q115 ), (|116|) , 
( |117| ) we see that ( |114j ) holds, and the proof is complete. □ 

Example: Assume that ab = p/q, 1 < p, q G Z, gcdfjo, q) = 1. If g is in the 
Feichtinger-algebra S [O, then by Theorem |372| /z° e S^. 



4 A Newton method for the computation of 
tight Gabor frames 

For applications and design purposes it is useful to have a fast algorithm 
for computing the canonical tight window h°. It is clear that in the case 
of integer oversampling this can be done via the Zak transform. However 
for non-integer oversampling the frame operator can at best only be block- 
diagonalized. This means that at some point we would have to compute 
the inverse square root of a matrix, which is computationally much more 
expensive than the inversion of a matrix. In this section we present a Newton- 
type algorithm that allows a very efficient calculation of h°. 

Newton's method is well-known in the context of computing matrix-sign 
functions and square roots of finite-dimensional matrices Jl8|, 24j. Let M be 



annxn matrix with positive real- valued eigenvalues and polar decomposition 
M = UH. Set M = M, X = M, and define the iterations 

Mfc+i^M. + tM*]- 1 !/), k = 0,l,... (119) 

X k+1 = \{X k + [X*]- 1 ), k = 0, 1, ... . (120) 

Then M k converges to and X k converges to U, see |18|, 23 1. 

In Section ^| we have seen that h° minimizes \\g — h\\ among all windows 
h for which (h, a, b) is a normalized tight frame. Clearly in the same way h° 
minimizes H7 — h\\. Also recall that 7 minimizes 

9 7 



:i2D 

\9\\ 1171 

among all dual windows 7. This suggests to compute h° iteratively by setting 
go = g and defining the iteration 

9w = i(M + M)' k = - 1 -- (122) 
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where 7° is the dual window associated with the frame (gk,a,b). Accord- 
ingly, we denote the frame bounds of the system (g^, a, b) by and Bk. 



Iteration ( 120 ) can be interpreted as a scaled Newton iteration. The proof 



of convergence of Newton's method as described in (p. 19 ) and ( |120|) relies on 



finite-dimensional methods (see 0]), which cannot be extended straightfor- 



wardly to infinite matrices or operators. 



4.1 Convergence of the algorithm 



In this section we analyze the convergence behavior of the iteration (|122|) . 
We need the following 

Lemma 4.1 Assume that (g, a, b) is a frame, and let (p be a function analytic 
in an open neighbourhood of [A, B], where A and B are the optimal frame 
bounds for (g, a, b), and assume that <p(s) > 0, s G [A, B] . Then (ip(S)g, a, b) 
is a frame with frame bounds 

min s(p 2 (s), max :s<p 2 (s). (123) 

A<s<B A<s<B 

Furthermore (g, a, b) and ((p(S)g, a, b) have the same canonically associated 
tight frame generating window, viz. h° . 

Proof: We have for any / G L 2 (M) that 

{f{S)g) nm )^{S)g) nm = p(S)[J2((v*(S)f,g nm )] = <p(S)S<p'{S)f. 

n,m n,m 

(124) 

Here we have used that 5*, and hence f(S), commutes with all relevant time- 
frequency shifts. Thus the frame operator corresponding to (ip(S)g, a, b) is 
given by ip(S)S<p*(S). By the spectral mapping theorem we have that 

<r[<p(S)S<p*(S)] = W(s)sif*(s) = s\ip(s)\ 2 \s G a(S)}, (125) 

and this gives the first part of the result. 

Next we compute the tight frame generating window h v canonically as- 
sociated to ((p(S)g, a, b) according to 

W = ^(S)S V (S)Y~ 2 v (S)g = S-H V 2 (S))- l ^(S)g = S^g = h°. (126) 

Here we have used that (p(s) > 0, s G [A, B] so that ip*(S) = tp(S) and 

_i 

<f(S) > so that ((</?(S)) 2 ) " = <p(S). □ 
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Theorem 4.2 Let (g,a,b) be a frame for L 2 (R) with optimal frame bounds 
A, B and canonical tight window h° . Set go = g and define 

ft=^ + rK). fc = l,2,.... (127) 



2Mbfc-ill IItLi 

i/ien g'fc converges quadratically to -^h°. 
Proof: Set 

gi = <p(S)g : p(x) = a + /fcr\ (128) 
where for simplicity we have denoted = ip-n- , j3k = jAm and a = ao, /? = A)- 

llflfcll llTfc II 

Note that 

A < — = J} 9 } - < B. (129) 
a \\S i g\\ 



Using Lemma |4.1| we obtain that 

Ai= min (a 2 x + 2a/3 + /3 2 x~ 1 ), (130) 

fix = max (« 2 x + 2a/3 + /5 2 x- 1 ). (131) 

xe[A,_B] 

As to Ai, we note that a 2 x + 2a (3 + /5 2 x _1 assumes its minimum value 4af3 
at x = (3 /a G [-4, -B], see ( |129| ). Hence Ai = 4a/3, and we furthermore note 
that 

4a ^MP^fOT^' (132) 

see (|39D. Hence Ai > -4, no matter what g is. As to B\ we note that 
a 2 x + 2af3 + /3 2 x _1 is strictly convex on [A, B] and assumes its minimum at 
x = (3 /a. Hence 

B x = max (a 2 x + 2af3 + f3 2 x~ 1 ). (133) 

x=A or B 

Therefore we find that 



Ax r 4a/3 4a/? 

mm 



B x Va 2 A + 2al3 + I3 2 A- V a 2 B + 2a(3 + ^B- 1 
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Figure 1: 



= min\ , 2;— L 134 

The graph in Figure Q shows the two functions /i(y) = (A 4 ^ )2 , / 2 (y) = 
for v = -. It is obvious that 

4AB <^< . (135) 



(A + E) 2 " B l ~ (y/A + y/B) 2 ' 



The above argument applies to any gk-i as in ( |127| ). In particular, we have 
for k = 1, 2, . . . 



1 " 7T " M + /? 15 = 77-^7 (136) 



Sfe_i 

Consider the recursion 



Ct = jrf^- * = 1 - 2 -- (137) 
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with C G (0, 1). The recursion ( |137| ) converges quadratically to 1, since the 
function 

f(y) = ^y 2 , ye (0,1] (138) 

increases on (0,1], satisfies /(0) = 0, /(l) = 1, /(y) > y,y G (0,1), and 
/'(l) = 0, /"(I) = —\- Evidently when we take Cq = A/B, we obtain that 

!>^>C k , k = l,2,... (139) 

and A k /B k — > 1, quadratically. 

It follows from (|127|) and (|1 28|) that the are all of the form (p k (S)g. 
Hence by Lemma |4.1| they all have h° as their canonically associated tight 
frame generator. At the same time, denoting by S k the frame operator of 

the system (g k ,a,b), we have that this canonically associated tight frame 

_ i 

generator is given by S k 2 g k . Hence 

h° = S;*g k , k = 1,2, (140) 

Since A k /B k — > 1, we easily see from 

g k = -(^+ S ^ 9k ~ 1 ) (141) 

that \\g k \\ — ► 1 as k— > oo. Also \\h°\\ = vab by (p9|), and it then follows 
from (|140[) and the fact that A k /B k — > 1 that 

A fc>J B fc ->— , A;^oo. (142) 



All these convergences are quadratic. Consequently S k 2 — » yabl quadrati- 
cally, which completes the proof. □ 

Considering ( |134j ) and the left-hand-side of (|135|) it is not difficult to 
check that quadratic convergence of the algorithm holds for other choices of 
the scaling parameters. Writing the iteration as 

9k+i = 2 ( a kg k + f3kli) , k = 0, 1, . . . (143) 
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it is easy to see that the algorithm converges quadratically to a scaled version 
of h° if ™ G [A, B]. The right-hand side in ( |135|) shows that the convergence 
cannot be expected to be faster than quadratic. Indeed, although the func- 
tion 

k(y)= n 4y 'x =/(y*), (144) 

(see (|138|) ) is considerably flatter than f(y) at y = 1, we have that fc(l) = 
l,fc'(l) = 0,fc"(l) = -§, hence we have at best that 1 - j*- = ^(1 - 

^-i) 2 compared to 1 - j*- = ±(1 - ^) 2 for /(y) in ( p8|) . This improves 
convergence behavior with respect to the constant, but not with respect to 
the order. However, since from a numerical viewpoint improving the constant 
is of considerable importance, we want to find out which choice of a&, leads 
to the optimal constant. 

Optimizing (|134|) with respect to - yields that the optimum choice for 
the scaling parameters is 

a k = (3~ l = {A k B k )~-^ (145) 

see also the graph in Figure |I|. The same optimal scaling parameter arises in 
the finite-dimensional Newton method, cf. JT7| as well in Balan's construction 
of "nearest tight frames", see Section 3 in 0. It is not difficult to see that 
for optimal scaling we get 

hence ^ takes on the maximum value in the right-hand side of (|135|) . 

The choice (|145|) however is practically not feasible, since it requires the 
computation (or estimation) of the frame bounds A k , B k in each iteration 
step, which is in general computationally expensive. 

An alternative is to find the scaling parameter a that minimizes 

Ho^-orVll- (147) 
A simple calculation shows that in this case 

IIt°II 

a = ill, (148) 

IMI 
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which is similar to Frobenius-norm scaling in (hence we will refer to it 
as Frobenius-norm scaling) and closely related to the scaling used in Theo- 



rem 



Next we estimate how the frame bounds get "tighter" from iteration k to 
k + 1 when using the scaling proposed in Theorem |4~2 . 



Corollary 4.3 Under the assumptions of Theorem [^.^ the frame bounds 
A k , B k of (g k ,a,b) can be estimated recursively by 



A k > 



1 2 V / Afc_i-Bfc_i 
ab v4 fc _i + B k _i 



B k < 



1 (Afe-i + flfc-i) 
ab 4A fc _i5 fe _i 



(149) 



Proof: It is sufficient to consider the case k = 1. We begin with the 
observation that for S satisfying AI < S < BI where A > 0, B < oo and 



there holds 



IS- 1 ! 



(150) 



Inequality ( |1 50| ) is equivalent to the operator version of the Kantorovich 
inequality ]2"2fl. For the reader's convenience we include a short proof. In- 



equality ( 150 ) follows from 



(A + B) 2 (S~ l x, x) 2 - AAB\\x\\ WS^xf 

1 2 



2AB\\S~ 1 x 



(A + B)\\S~*x\ 



+ 



+4AB((BI -S)(I- AJS-^S-tx, S~*x) WS^xW > 0, 

since (BI — S)(I — AS' 1 ) is a non-negative operator. 
As a consequence we have now immediately that 



(151) 



At 



1 



(S- l g,g) > 1 2VAB 



\g\\\\S 1 g\\ ab\\g\\\\S 1 g\\ abA + B 



(152) 



which is the first inequality in (|149 ). 
From ( |135| ) we have 



1 (A + B) 2 (S- l g,g) ^ 1 (A + B 



AAB ab AAB WgUS^g] 

and this is the second inequality in (HU 



> 



ab AAB 



(153) 

a 
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Both estimates in ( |149j ) are sharp. This can be seen by considering the 
following example. Assume a = b = 1. Let e > be arbitrarily small, and 
take g such that \Zg\ = B^ ,l/\Zg\ = \ZS~ l g\ = B~? on a set of measure 
1 — e in the unit square while \Zg\ = , \j\Zg\ = \ZS~ x g\ = A~% on the 
complement of this set. Then we have \\g\\ ~ , HS 1-1 ^ « B~z, whence 

B l = max + - \ : .. + ttt| : (154) 

/x 1 or 1 \ (A + 5) 2 . . 

k max 1 1 =- — '—. (155 

x=a,b\AB 2 AB- 1 ) AAB v ; 

Hence the second inequality in ( |149| ) is sharp. 

As to sharpness in the first inequality in (|149Q , we consider a g that has 
a Zak transform with Zg = B? on a set C [0, l) 2 of measure t G [0, 1] and 
Zg = A? on the complementary set. Then 

Al = c-l II = My ||||tc-l || ( 156 ) 

\\9\\\\S 1 g\\ \\Zg\\\\ZS 1 g\\ 
(A + t(B-A))(j + t(±-j))y l2 (157) 

Taking t — \ here, we have 

V AAB ) A + B' K J 

which is the right-hand side of the first inequality in ( 149|) . 
We compute Zg\ for this last example. We find 

i i 

|( a+4-A) ) 2 + l( ^(B-A) ) 2 whenZ(7 = Ai 



The difference of upper and lower number in ( |160| ), ignoring terms of order 
(2t — l) n (B — A) m with n + m > 4 (there are no terms with n + m = 4), is 

2VAB 
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E.g., when t = | the difference is zero. This example also illustrates that the 
algorithm may converge even much better than predicted by the theoretical 
estimates presented above. 



4.2 Numerical simulations 

The numerical simulations presented in this section are representative for 
a large number of situations. We consider a well-localized window g and 
oversampling rate -4 G [1,2]. Such a setup is typical for oversampled modu- 
lated filter banks as well as for OFDM systems (where fc-times oversampling 
corresponds to fc-times undersampling) . We choose g(t) = 2^e~ n ^ 2 with 

a = b = y^jf resulting in an oversampling rate of 16/15. For the numerical 

implementation we follow the finite-dimensional model in [Q, Sec. 1.6. 

In the first experiment we consider different types of scaling for itera- 



tion ( |143| ). We compare the Newton method with norm-scaling as used in 



Theorem |4.2j to optimal scaling (cf. ( |145|) ), Frobenius-norm scaling (cf. Q148Q ) 



and to the unsealed Newton iteration (i.e., = (3k — 1). We precompute 
h°, calculate in each iteration the normalized error ||<7fc/||<7fc|| — ^ /||^°|||| an d 
terminate the iterations when the accuracy is within 10~ 14 . The results are 
shown in Figure 0. 

As expected, optimal scaling results in fastest convergence. Norm-scaling 
and Frobenius-norm scaling yield almost identical convergence, slightly slower 
than optimal scaling. Recall however, that the costs for computing the opti- 
mal scaling parameters are much larger than one iteration step for the other 
two scaling schemes. Thus the faster convergence comes at a prohibitive high 
price. Newton's method without scaling requires significantly more iterations 
in order to achieve the same approximation error. In terms of overall compu- 
tational costs Newton's method with norm-scaling or Frobenius-norm scaling 
are clearly the most efficient methods. 



An alternative to the algorithm presented in Theorem (12 may be to 
compute first iteratively and at the end calculate h° = S~^g. Efficient 
iterative methods to compute the (inverse) square root of a positive definite 
matrix have been proposed in P5fl . These methods are different from 
Newton's method as described in (|119| ). 



Let M be a symmetric positive definite matrix. Then Sherif 's method to 
compute M~2 is given by (cf. p0|): 

Mn = I 
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Figure 2: Comparison of different scalings for Newton's method. 



M k+l = 2M k (I + MM 2 k )-\ fc = 0,l,... (162) 

This method yields quadratic convergence. The method due to Lakic achieves 
even a cubic rate of convergence. Its iteration rule is (cf. |25|): 

M = I 

M k+1 = l -M k {I + 8(1 + 3MA4 2 )- 1 ), k = 0, 1, . . . (163) 

Lakic's method yields usually slightly faster convergence than our method, 
whereas Sherif's method is somewhat slower than ours, see Figure [|. For 
an overall comparison of the efficiency of iterative methods we have to take 
into account the computational effort needed for one iteration step. The 
number of flops for the matrix-multiplications and inversions in ( |141| ), ( |162| ), 



and ( |163|) depends on the chosen representation of the frame operator. Note 
that for the proposed method this representation has to be computed for each 
iteration (since we are dealing with a new frame in each iteration), which is 
not the case for the other two methods. However Newton's method needs only 
one matrix inversion per iteration, whereas the other two methods require 
three matrix-matrix multiplications in addition to one matrix inversion. 
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Newton w. norm-scaling 
Sherif's method 
Lakic's method 




Figure 3: Comparison of different methods to compute the tight window h c 



Fortunately the computation of the various representations of the frame 
operator can in general be done very efficiently For the representation in 
the time-domain this amounts to simple data addressing, in the other three 
domains FFT-based algorithm can be used (possibly combined with data ad- 
dressing). Thus, as long as the computation of this representation is cheaper 
than the inversion of the resulting matrix, the method proposed in Theo- 



rem [4.2| outperforms the other two methods in terms of computational effi- 
ciency 
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